全 极 化 GNSS-R 陆 面 参数 延迟 多 普 勒 图 模型 
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Abstract: GNSS-R is a new promising remote sensing technique, which utilizes the reflected signals of GNSS 
constellation to remotely sense the ocean or land surface. At present, most of the existing research are based on 
the experimental observations, few theoretical studies have been carried out. This paper has developed a fully 
polarization GNSS-R delay Doppler map for land geophysical parameters study. This theoretical model was based 
on the ocean surface GPS scattering model, after the modifications, it has been used for the land surface. As for 
the calculations of land geophysical parameters (bare soil and vegetation), random surface scattering model and 
the first order radiation transfer equation model were used. Effects of bare soil and vegetation parameters on 
the delay Doppler maps were simulated. In order to use the polarization information for backward inversion, our 
developed theoretical model has the ability of fully polarizations calculations, i.e. RR, LR, HR and VR pol. This 
theoretical model is a mechanism tool for the data explanation of space borne mission, experimental campaign 


design, data simulations and the backward inversion algorithm development. 
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摘要 : 
目前 ， 针 对 GNSS-R 的 研究 多 从 观测 角度 出 发 ， 对 其 散射 机 理 研 究 相 对 较 少 。 本 文 针 对 GNSS-R 


(GNSS-Reflectometry) 反射 信号 陆 面 参数 的 延迟 多 普 勒 图 理论 模型 展开 研究 ， 对 于 陆 面 参数 〈 裸 士 、 植 
WO. 的 双 站 散射 特性 的 计算 采用 修改 后 的 随机 粗糙 面 散射 模型 和 零 阶 / 一 阶 辐射 传输 方程 模型 。 通 过 该 模 
型 模拟 分 析 了 裸 士 、 植 被 参数 变化 时 ， 相 应 的 延迟 多 普 勒 图 的 变化 。 同 时 为 利用 极 化 信息 进行 相应 地 物 
参数 的 后 向 反 演 ， 模 型 中 增加 全 极 化 计算 功能 。 全 极 化 GNSS+R 陆 面 参数 延迟 多 普 勒 图 理论 模型 的 研究 有 
助 于 星 载 数据 的 解释 、 地 表 实 验 设计 、 数 据 仿真 和 后 向 反 演算 法 的 开发 。 
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星 的 反射 信号 对 地 物 参数 遥感 的 新 兴 对 地 观测 方式 ， 在 国内 外 引起 广泛 关注 
[1.2]。 由 于 海洋 表面 相对 均一 ， 极 化 特性 不 明显 ， 因 此 GNSS-R 在 海洋 领域 的 研究 相对 较 早 ， 较 成 熟 
[4,5]。 针 对 陆 面 ， 目 前 国内 外 研究 主要 依赖 于 地 面试 验 手 段 ， 建 立 GNSS-R 信号 或 GPS 多 路 径 数 据 与 地 表 
参数 间 的 相关 性 ， 以 研究 其 应 用 的 可 行 性 或 者 建立 区 域 性 的 定量 反 演算 法 。 
早期 的 SMEX02 实验 中 首次 包含 机 载 GPS-R 土壤 水 分 实验 ， 接 收 机 天 线 仍 然 指向 天 底 ， 但 增加 了 接收 
机 天 线 增 益 ， 以 期 获得 更 好 的 SNR 数据 ， 该 实验 验证 了 GPS 反射 信号 和 土壤 介 电 常数 以 及 土壤 体积 含水 
量 之 间 的 关系 ， 利 用 GPS 数据 获取 的 介 电 常 数 和 Wang-Schmugge 理论 模型 吻合 较 好 ， 但 是 均 较 低 于 实测 
土壤 水 分 数据 [5]。 

随后 的 BAO 塔 实验 [6]， 接 收 机 天 线 包 含 多 种 天 线 类 型 . 低 增益 半球 LHCP 天 线 ;，4 个 高 增益 (~12dB) 
V、H、LHCP (Right Hand Circular Polarization) 和 RHCP 天 线 ; 为 了 分 析 极 化 特性 对 反射 信号 反 演 的 影响 ， 
采用 一 阶 SSA (Small Slope Approximation) 模型 计算 土壤 散射 特性 ， 但 是 由 于 过 于 简单 的 模型 假设 造成 
BAO 塔 实验 中 的 实测 数据 无 法 验证 相应 的 理论 模型 。 

GPS-IR 方法 不 需要 研制 专门 接收 机 ， 直 接 采 用 地 球 物理 或 者 地 球 测绘 类 接收 机 ， 即 利用 其 多 路 径 信息 
(直射 信号 和 反射 信号 的 相干 信号 ) 对 地 物 参数 进行 遥感 监测 [7]。 采 用 电磁 单 次 散射 前 向 模型 验证 近 地 
表土 壤 水 分 和 SNR 相关 量 的 变化 关系 [8]， 研 究 表明 地 表 5cm 以 上 区 域 的 土壤 水 分 变化 会 影响 有 效 反射 计 
高 度 、 相 位 和 幅度 ， 浅 层 (<1-cm DO 的 土壤 水 分 变化 对 其 影响 最 大 ， 土 壤 质 地 类 型 的 影响 可 以 忽略 ， 相 
位 是 土壤 水 分 变化 反 演 的 最 有 效 度量 标准 [9]。 

IPT 方法 是 利用 测 得 的 直射 信号 和 反射 信号 的 相干 能 量 反 演 地 形 、 植 被 覆盖 区 土壤 水 分 和 植被 高 度 参 
数 [10]， 该 种 方法 中 接收 到 的 能 量 信号 中 ， 四 槽 位置、 振幅 以 及 个 数 与 地 物 参数 有 关 [1112]， 实 际 应 用 中 
利用 其 信息 进行 反 演 : 地 形 反 演 精度 是 RMSE=27cm; 植被 覆盖 区 土壤 水 分 的 精度 为 2%~5%; 植被 高 度 的 
反 演 精度 为 RMSE=3~5cm。 


GNSS-R 的 研究 目的 是 实现 星 载 观测 ,美国 NASA 的 8 星星 座 “ 飓 风 全 球 导 航 卫 星系 统 ” CYGNSS 


GNSS-R 是 利用 导航 


(Cyclone Global Navigation Satellite System) 已 经 于 2016 年 12 月 12 日 成 功 发 射 ， CYGNSS 星座 发 射 的 主 
要 目的 是 利用 GPS 直射 信号 和 反射 信号 测量 热带 风暴 和 朵 风 期 间 海 洋 表 面 风 场 信 息 ， 同 时 CYGNSS 也 为 陆 
面 参数 及 冰雪 圈 研 究 提 供 了 重要 契机 [13]， 其 成 功 发 射 标志 着 GNSS-R 星 载 计划 走向 业务 化 ，CYGNSS 的 


接收 机 为 DMR (Delay Mapping Receiver)， 每 个 观测 平台 中 包含 一 个 DDMI (Delay Doppler Mapping 
Instrument) 接收 仪 ， 接 收 机 获得 的 是 海面 或 陆 面 参数 的 DDM 波形 ， 而 DDM 波形 是 GNSS-R 技术 反 演 地 


球 物理 参数 的 基本 观测 量 ， 其 仿真 软件 是 开展 GNSS-R 基本 理论 研究 和 项 目 工程 参数 设计 的 重要 工具 。 


针对 星 载 观测 的 仿真 软件 如 CYGNSS E2ES 软件 [14] ( Brien et al., 2017), P2EPS 仿真 软件 [15](Park et al., 
2017), GREEPS 仿真 软件 [16](Bai et al., 2016) 均 针对 海洋 表面 参数 建立 。 目 前 针对 对 陆地 表面 的 DDM 仿真 
系统 较 少 。 


本 文 目 的 是 发 展 建立 针对 不 同 地 物 参 数 〈 裸 土 、 植 被 ) 的 全 极 化 GNSS 反射 信号 模型 。GNSS 直射 信 
号 经 地 表 反 射 后 ， 被 专门 的 机 载 / 星 载 GNSS-R 接收 机 所 接收 ， 形 成 一 种 双 站 /多 站 雷达 的 工作 模式 。 现 有 
的 地 物 参数 微波 散射 模型 ， 多 是 针对 单 站 雷达 的 后 向 散射 模型 或 者 是 针对 辐射 计 的 被 动 微波 反射 率 模型 
[17]。 因 此 针对 GNSS-R 遥感 ， 需 要 对 现 有 常规 微波 散射 模型 进行 修改 ， 使 得 模型 具备 双 站 散射 计算 能 力 
[18]。 

与 传统 单 站 雷达 或 者 辐射 计 不 同 ， 导 航 卫 星 发 射 的 为 圆 极 化 信号 (如 GPS 卫星 发 射 的 为 RHCP)， 该 
信号 经 地 表 反 射 后 ， 极 性 发 生 改 变 。 由 于 电磁 波 的 极 化 特性 受 目 标 物 介 电 常数 、 物 理 特性 、 几 何 形状 和 
取向 等 因素 影响 ， 对 GNSS-R 极 化 特性 的 研究 有 利于 探测 目标 参数 提取 ， 因 此 本 文 发 展 建立 了 全 极 化 散射 
模型 。 

在 上 述 发 展 建立 的 双 站 全 极 化 模型 的 基础 上 ， 发 展 建立 了 针对 典型 陆 面 参数 土壤 水 分 、 植 被 ) 的 

- 全 极 化 DDM 仿真 模型 。 

rs 文章 中 第 二 部 分 对 发 展 建立 的 全 极 化 GNSS-R 反射 信号 陆 面 参数 延迟 多 普 勒 图 模型 的 具体 理论 和 研究 
方法 进行 介绍 ， 包 括 双 站 全 极 化 散射 模型 和 陆 面 GNSS 反射 信号 模型 ， 即 拓展 的 Z-V 模型 以 及 相应 的 工作 
流程 ;第 三 部 分 则 给 出 了 相应 的 模拟 结果 ， 包 括 裸 土 、 植 被 相关 参数 变化 时 ， 相 应 的 各 种 极 化 双 站 雷达 
散射 截面 的 变化 关系 以 及 对 应 的 延迟 多 普 勒 图 的 变化 关系 。 第 四 部 分 对 结果 进行 讨论 ， 第 五 部 分 为 总 


Ql ”2、 理 论 及 方法 
> 本 部 分 对 全 极 化 GNSS-R 散 射 信号 陆 面 参 数 物理 模型 的 理论 基础 和 工作 流程 进行 介绍 ,包括 计算 裸 土 、 
ci o0 植被 覆盖 地 表 的 双 站 全 极 化 散射 模型 ， 计 算 相应 地 物 参数 的 GNSS-R 散射 信号 延迟 多 普 勒 图 模型 以 及 该 理 
Io 论 模型 的 总 体 工作 流程 ， 
O — 2a 双 站 全 极 化 散射 模型 

对 于 裸 土 散射 特性 的 计算 分 为 两 部 分 ， 假 设 地 表 足 够 光滑 ， 可 以 用 光滑 表面 的 反射 率 模型 对 地 表 的 

散射 特性 进行 描述 。 对 于 非 相干 部 分 的 计算 ， 采 用 随机 粗糙 面 散射 模型 [19]。 
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qp 


HH, O° 表示 双 站 雷达 散射 截面 ， 下 标 p. q 分 别 代表 发 射 和 接收 的 极 化 状态 。 该 模型 整体 上 分 为 基 尔 


霍 夫 项 ce ， 基 尔 霍 夫 项 的 补偿 项 CQ 名和 两 者 的 交叉 项 Cg 三 项 之 和 (公式 1)。 


对 于 植被 散射 特性 的 计算 采用 迭代 算法 求解 辐射 传输 方程 [17]: 假定 向 上 和 向 下 传输 的 能 量 强度 矩阵 
分 别 为 "，[r， 可 以 通过 辐射 传输 方程 的 一 阶 微分 形式 表示 ， 如 公式 2-3. 


A nps eri (2) 
dz : 
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—I (z)-k,r (z)-F (z) (3) 
dz 
式 中 
k„ =k, /cos0, = k,/ 4, (4) 
44, =cos0,, u =cos0 (5) 
I*(z)=1(z,4.) (6) 
I (z)-I(z.-u.d.) (7) 


式 中 对 角 和 矩阵 大 k, 分 别 为 体 散 射 系 数 矩 阵 和 消光 系数 年 阵 ， 大 = 大 /cosg , Pl 14,468 0) wit 


射 相 位 矩阵 ， 辐 射 传输 方程 解 的 一 阶 积分 方程 形式 如 下 : 
r'(z-r (-aje** «p (zee le (8) 


r (zer (oye +f F (z e ae (9) 
ser F^ 为 向 上 向 下 传输 的 能 量 强度 的 函数 。 
$ B k 2m 
F(z)-. 


k 2z el ] 
e [,P(5u..74.6, -0)1 d ud 


f Peu. usd. - à) dude 


0 
(10) 


针对 植被 的 辐射 传输 方程 模型 [20,21], 假定 入 射 强度 为 7(6,g) ,散射 强度 为 1 (8...) 其 中 加 ,外 分 别 


是 入 射 的 天 项 角 和 方位 角 ， ,gp 为 散射 的 天 顶 角 和 方位 角 。 
I, (8..9.)- —L., (0,,0,:0,,0,50,,0,) 1, (6. 9,) S 
入 射 能 量 被 单 散 射 体 散射 后 , SURE Be tech FE E NET Be tec os BE Z [8] n] DA SEE PRU] 2 UAR E Lm 联系 ( 公 
X 1D. (0,,0,) 为 单 散射 体 朝向 ，r 为 散射 强度 与 散射 体 之 间 的 距离 ，Lnm AI AS A FL NT ABE, HE 
射 振幅 S 定义 。 


c, (v, x uz) 40¥ MY; (12) 


对 于 各 种 极 化 散射 截面 的 计算 采用 极 化 合成 的 方法 [22] (公式 12)， 式 中 三 ! ALY” 分别 为 发 射 机 和 接 


收 机 的 归 一 化 斯 托 克 斯 矢量 ， 通 过 对 其 中 椭 倾 角 和 椭 率 角 的 变换 ， 可 以 得 到 任意 极 化 的 雷达 双 站 散射 截 
面 。 


2.2 全 极 化 GNSS-R 散射 信号 模型 
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早期 的 z-v 模型 是 应 用 海洋 表面 的 相关 函数 和 相关 功率 模型 [23]， 本 质 上 是 双 站 雷达 的 积分 形式 ， 是 
昌 收 发 分 置 的 散射 问题 近似 派生 出 来 的 , 该 模型 最 初 是 用 来 研究 海洋 表面 风 场 , 由 基 尔 霍 夫 光学 模型 出 发 。 
本 文 发 展 建立 的 理论 模型 将 其 拓展 用 于 裸 土 、 植 被 延迟 多 普 勒 图 仿真 分 析 。 重 点 对 其 积分 区 域 、 观 测 几何 
以 及 观测 目标 散射 模型 进行 修改 [24]。 该 模型 可 以 描述 在 不 同时 间 延 迟 和 多 普 勒 频 移 下 GPS 反射 信号 的 能 
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Y, (2,7)= TPA I d (£ - z)sinc? (/ - f)a4 (13) 


式 中 起， 接收 到 的 GPS RAMSAR, Rae? 和 频率 了 Home, Pn. J GPs 卫星 发 射 的 能 量 ; 
Gr: 发 射 机 天 线 增 益 ，Cx: PML ABM A; Res 接收 机 到 地 面 反射 点 的 距离 ， 心 : 发 射 机 到 地 表 镜 像 
RAWE: A: 电磁 波 波 长 ;7 : 信号 处 理 时 采用 的 相干 积分 时 间 ; o^: 无 维 ， 双 站 雷达 散射 截 


Wi A(f-r): GPS 相关 函数 (三 角 函数 ) , 7,7 分别 为 复制 信号 和 入 射 信号 延迟 ，sinc*[ 了 一 :多 普 


勒 滤波 函数 ，f, f 分别 为 复制 信号 和 入 射 信号 的 频率 ，A: 有 效 散射 面积 ， 接 近 闪 烁 区 ，dA:A 区 域内 的 
积分 面积 。 
2.3 理论 方法 流程 图 

GNSS-R 反射 信号 地 物 参数 理论 模型 工作 框架 如 图 1 所 示 ， 其 中 裸 土 和 植被 的 双 站 雷达 散射 截面 分 别 通 
过 随机 粗糙 面 散射 模型 和 一 阶 辐射 传输 方程 模型 计算 得 到 ， 由 GPS 卫星 及 相应 的 GNSS-R 接收 机 的 位 置 和 
速度 信息 计算 得 到 卫星 几何 模型 中 所 需 参 数 ， 最 后 由 全 极 化 GNSS 散射 信号 模型 计算 得 到 相应 地 物 参数 的 
延迟 多 普 勒 图 。 


全 极 化 GNSS 散射 信号 模型 
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图 1. 全 极 化 GNSS+R 陆 面 参数 延迟 多 普 勒 图 理论 模型 工作 流程 图 
Fig.1 Flowchart of the fully polarimetric GNSS+R delay Doppler map model for land geophysical parameters 


3 模拟 结果 
在 以 下 模拟 中 GPS 载波 频率 采用 L1 波段 ， 积 分 时 间 为 0.001s。GPS 卫星 (tx) 和 相应 GNSS-R 接收 机 


(rx) 位 置 矢量 (_pos 单位 km) 和 速度 矢量 (_vel 单位 m/s) 信 息 如 下 : 


tx pos = [-11179 -13160 20342]; 
tx vel - [2523 -361 1163]; 
rx pos = [-4070 -3583 4528]; 
rx vel = [-4738 -1796 -5655]; 
3.1 裸 土 、 植 被 DDM 波形 
地 表 模 型 输入 参数 〈 地 表 和 植被 ) 如 表 1 所 示 。 通 过 输入 不 同 的 植被 和 土壤 参数 可 以 获取 其 对 双 站 
雷达 散射 截面 和 最 终 的 DDM 波 
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# 1. GNSS-R DDM 地 表 模 型 输入 参数 。 树 冠 层 高 度 为 11.00 米 ， 冠 层 密度 为 0.2 R/M 。 
Tab.1 Model inputs for GNSS-R DDM simulator. Canopy height is 11.00 m, canopy density is 0.2 trees/m" 


树冠 层 

树干 枝条 叶子 EPIRI 
介 电 常数 (D 16.4+7.3]  16473j 19+6j 27-12.43j 
重力 含水 量 (g/cm) 0.6 0.6 0.8 0.8 
于 物质 密度 (gem?) 0.4 1.0 0.976 0.1 
长 度 (cm) 1600 200 0.02 1.7 
半径 (cm) 20.08 2.0 4.0 0.1 
密度 (散射 体 /m3) 0.2 3.4 85000 85000 
地 表层 
介 电 常数 〈/) 6.3+1.6j 
均 方 根 高 度 (cm) 0.45 
相关 长 度 (cm) 18.75 


由 三 角 函 数 人 和 多 普 勒 滤波 函数 3 分 别 确定 的 环形 区 和 多 普 勒 等 值 区 如 图 2 所 示 。 
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UR 1 中 相应 参数 作为 模型 的 输入 ， 在 


2 Illustration for time delay and Do 
图 2 所 示 的 时 间 延 迟 和 多 普 勒 频 移 下 ， 植 被 的 全 极 化 双 站 
散射 截面 示意 图 如 图 3 所 示 。 图 3 分 别 给 出 了 不 同 极 化 RR、LR、VR 和 HR) 的 双 站 雷达 散射 截 


图 2， 时 间 延 迟 和 多 普 勒 频 移 示意 图 
ppler shift. 


图 可 以 看 出 ， 在 中 心 区 域 ， 即 镜像 散射 点 附近 双 站 雷达 散射 面 截 较 强 ， 而 在 


图 。 从 双 站 雷达 散射 截面 


漫 散射 区 域 ， 相 应 的 雷达 散射 特性 较 弱 。 相 同时 延 和 多 普 勒 频 处 ， 不 同 极 化 的 散射 波形 差异 较 大 ， 因 


为 在 不 同 极 化 时 ， 地 物 参数 的 散射 特性 差异 较 大 。 


Fig.3 Vegetation bistatic radar cross section fo 


当 模 型 中 去 掉 树 冠 层 ， 只 保留 地 表层 ， 即 裸露 地 表 的 双 站 雷达 


射 模型 计算 得 到 。 图 4 | 
散射 截面 差别 较 大 。 由 于 


图 3 植被 RR、LR、VR 和 HR 极 化 双 站 雷达 散射 截面 图 


RR. LR, VR and HR pol. 


散射 截面 的 计算 采用 随机 粗糙 面 散 


给 出 植被 和 裸 土 全 极 化 双 站 雷达 散射 截 


散射 点 对 应 的 观测 几何 的 差异 ， 


应 的 双 站 雷达 散射 截面 差 不 同 ， 由 于 植被 的 衰减 ， 植 被 双 站 雷达 散射 截面 小 于 裸 土 的 双 站 雷达 散射 截 


用 ， 具 体 衰减 程度 与 相应 


的 观测 角度 有 关 。 


差别 示意 图 ， 不 同 极 化 的 双 站 雷达 


Ear 


在 各 种 极 化 的 不 同 观测 几何 下 ， 相 
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Fig.4 Bistatic radar cross section difference map for vegetation and bare soil 


图 5 则 是 在 不 同 极 化 时 ， 对 应 的 植被 和 裸 土 延 迟 多 普 勒 散射 信号 能 量 差 示意 图 。 不 同 极 化 下 ， 散 


射 信号 能 量 差 别 较 大 。 因 此 可 以 有 效 利用 极 化 信息 进行 地 表 和 植被 相关 参数 的 监测 和 反 演 。 由 于 地 物 
参数 的 散射 特性 受 观测 几何 ( 即 散 射 的 天 顶 角 和 方位 角 〉 影响 较 大 ， 因 此 散射 信号 能 量 差 在 不 同 观测 


几何 下 呈现 不 同 特性 。 
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图 5. 植被 覆盖 与 裸 土 全 极 化 延迟 多 普 勒 图 差别 图 


Fig.5 Delay Doppler difference map for vegetation and bare soi 


针对 裸 土 ， 该 理论 模型 可 以 模拟 分 析 地 表土 壤 水 分 、 地 表 粗 糙 度 变化 对 延迟 多 普 勒 图 的 影响 ， 针 对 植 
被 ， 该 理论 模型 可 以 模拟 分 析 植 被 生长 状况 、 生 物 量 变化 、 植 被 含水 量 以 及 地 表 参 数 变化 对 延迟 多 图 普 
勒 图 的 影响 。 


4、 讨 论 
目前 的 微波 散射 模型 ， 多 是 针对 辐射 计 的 发 射 率 模型 或 者 是 针对 单 站 雷达 的 后 向 散射 模型 。GNSS-R 
遥感 中 GNSS 导航 卫星 群 和 相应 的 反射 信号 接收 机 之 间 构成 典型 的 双 站 雷达 工作 模式 ， 因 此 计算 时 需要 采 
用 双 站 雷达 模型 ， 即 模型 需要 可 以 计算 任意 散射 方向 (散射 天 顶 角 和 方位 角 ) 的 散射 模型 。 本 文 计算 神 
土 用 的 高 级 积分 方程 模型 和 计算 植被 的 辐射 传输 方程 模型 均 为 双 站 雷达 散射 模型 ， 在 进行 模型 验证 时 
将 该 模型 的 角度 〈 散 射 天 项 角 和 方位 角 ) 修改 为 后 向 单 站 雷达 的 角度 ， 得 到 的 模型 计算 结果 与 未 经 修改 
的 模型 计算 结果 相同 ， 以 此 表明 各 种 角度 〈 散 射 天 项 角 和 方位 角 ) 的 双 站 散射 计算 能 力 的 正确 性 。 
为 降低 电离 层 对 GNSS 卫星 信号 的 影响 ， 导航 卫星 群 作为 发 射 源 发 射 右 旋 圆 极 化 的 信号 ， 因 此 需要 
双 站 散射 模型 可 以 计算 包含 圆 极 化 的 全 极 化 的 计算 功能 ， 通 过 极 化 合成 公式 可 知 ， 归 一 化 斯 托 克 斯 矢量 
中 机 倾角 和 机 率 角 决定 着 不 同 极 化 方式 ， 针 对 全 极 化 的 计算 功能 ， 将 极 化 合成 公式 中 的 椭 倾 角 和 椭 率 角 
修改 为 线 极 化 的 角度 ， 将 发 展 建立 的 模型 计算 结果 与 原 有 未 经 修改 的 线 极 化 模型 的 计算 结果 进行 比较 ， 
二 以 此 验证 全 极 化 模型 计算 的 正确 性 。 
>< 本 文 发 展 建立 的 全 极 化 ODM 仿真 模型 是 在 原 有 海洋 表面 DOM 仿真 模型 的 基础 上 发 展 建立 的 ， 模 型 
验证 时 将 双 站 雷达 散射 截面 用 相应 的 土壤 水 分 或 植被 的 双 站 雷达 散射 截面 来 蔡 代 ， 模 拟 结果 的 正确 性 验 
Tao 证 时 ， 需 要 验证 原 有 海洋 表面 DDM 仿真 模型 的 双 站 雷达 散射 截面 接口 设置 是 否 合理 ， 因 此 将 海洋 表面 的 
双 站 雷达 散射 截面 作为 接口 处 输入 值 ， 计 算得 到 结果 与 原 有 模型 结果 相同 ， 以 此 验证 模型 的 正确 性 。 
作为 对 地 观测 新 方式 ，GNSS-R 遥感 的 发 展 在 国内 外 引起 广泛 关注 。GNSS+R 是 指 除了 利用 导航 卫星 的 反 
射 信号 以 外 ,还 可 以 利用 其 他 数字 通讯 卫星 的 反射 信号 对 地 物 参数 进行 遥感 。 GNSS-R 遥感 的 一 个 特点 就 是 
不 断 变 化 的 观测 几何 ， 而 且 其 工作 波段 只 有 上 波段 。 利 用 其 反射 计 的 思想 ， 采 用 数字 通信 卫星 的 反射 信号 
为 例 对 海面 粗糙 度 进行 研究 实验 上 验证 了 其 应 用 的 可 行 性 [25,26] 。 由 于 数字 通信 卫星 工作 在 微波 频段 ,所 
以 从 理论 上 进行 土壤 水 分 监测 和 植被 特性 研究 是 可 行 的 。 但 有 必要 对 其 散射 机 制 进行 研究 ， 明 确 机理 ， 并 
展开 相应 的 实验 。 
本 文 发 展 建立 的 双 站 全 极 化 散射 模型 频率 覆盖 范围 为 0.5-10GHz 观 测 入 射 天 项 角 大 于 10”， 在 针对 相 
应 的 导航 卫星 或 者 其 他 数字 通讯 卫星 的 三 角 函数 和 多 普 勒 滤波 函数 的 修改 基础 上， 该 模型 可 以 拓展 到 
GNSS+R 的 应 用 。 本 文中 散射 特性 的 研究 为 GNSS+R 相关 接收 机 设计 、 数 据 仿真 和 定量 反 演算 法 的 开发 可 以 


提供 关键 的 机 理工 具 , 而 给 出 的 陆 面 参数 ( 裸 土 、 植 被 、 冻 融 地 表 ) 的 延迟 多 普 勒 图 仿真 分 析 模 型 为 GNSS+R 
的 发 展 提供 了 一 定 的 理论 基础 。 


5、 结 论 

GNSS-R 工作 模式 本 质 上 为 双 站 雷达 ， 接 收 机 接收 到 的 数据 为 延迟 多 普 勒 图 。 本 文 发 展 建立 了 全 极 化 
GNSS-R 陆 面 参 数 延 迟 多 普 勒 图 理论 模型 。 该 理论 模型 基于 海洋 表面 的 GNSS 散射 信号 模型 ， 针 对 陆 面 的 
裸 土 、 植 被 应 用 展开 。 相 应 地 物 参数 双 站 雷达 散射 截面 的 计算 分 别 采 用 修改 后 的 随机 粗糙 面 散射 模型 、 
零 阶 /一 阶 辐射 传输 方程 模型 。 利 用 极 化 合成 的 方法 使 得 微波 散射 模型 可 以 计算 各 种 极 化 的 散射 特性 。 当 
地 表 有 无 植被 覆盖 以 及 地 表 参 数 变化 时 ， 会 导致 各 种 极 化 下 的 雷达 散射 截面 差异 和 延迟 多 普 勒 图 的 变 
化 。 该 理论 还 可 以 用 来 计算 土壤 水 分 、 植 被 生长 状况 ， 生 物 量 或 者 含水 量变 化 时 ， 延 迟 多 普 勒 图 的 变化 
情况 。 该 延迟 多 普 勒 图 理论 模型 可 以 计算 各 种 极 化 下 ， 相 应 地 物 参数 的 延迟 多 普 勒 图 的 变化 ， 有 利于 利 
用 极 化 信息 来 提取 相应 地 物 参数 的 有 效 信息 。 基 于 微波 散射 机 理 的 延迟 多 普 勒 模型 的 理论 研究 工作 有 助 
于 星 载 GNSS-R 观测 数据 的 解释 和 分 析 ， 并 有 利于 后 向 反 演 算法 的 开发 ， 其 研究 结果 也 为 GNSS+R 陆 面 遥 
感 的 发 展 莫 定 一 定 的 理论 基础 。 
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